The trumpet solution from spherical gravitational collapse with puncture gauges 
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We investigate the stationary end-state obtained by evolving a collapsing spherical star with the gauges rou- 
tinely adopted to study puncture black holes. We compare the end-state of the collapse with the trumpet solution 
found in the evolution of a single wormhole slice and show that the two solutions closely agree. We demonstrate 
that the agreement is caused by the use of the Gamma-driver shift condition, which allows the matter to fall in- 
wards into a region of spacetime that is not resolved by the numerical grid, and which simultaneously finds the 
stationary coordinates of the trumpet outside the matter. 



I. INTRODUCTION 

The possibility of obtaining stable evolutions of relativis- 
tic compact objects within the framework of numerical rel- 
ativity relies crucially on the choice of suitable coordinate 
gauge conditions. In particular, the well-known puncture ap- 
proach relies on the gauge conditions to deal with explicit 
coordinate singularities present in the computational domain. 
These gauge conditions are built upon the "1+log" lapse con- 
dition 0]] and the "Gamma-driver" shift condition JH, and 
allow the coordinate singularity (puncture) to be advected 
across the grid by the shift fl3H7]]. Hereafter we will refer to 
these gauge conditions as the puncture gauges. 

It has been shown SIShU that, for a single puncture with- 
out linear momentum or spin, the wormhole topology of the 
puncture initial data foil ceases to be resolved by the mesh 
in the code. Instead, the mesh approaches an asymptotically 
cylindrical stationary, solution, i.e., the trumpet solution. The 
puncture gauges can thus be viewed as a sort of "natural ex- 
cision" which squeezes the singularity into an unresolved re- 
gion of the computational domain. It is nowadays commonly 
used in binary black hole simulations by many groups follow- 
ing si. 

The puncture gauges have furthermore been tested in a 
number of codes solving also the equations of relativistic hy- 
drodynamics H4I - H9I1 . It has been found effective in han- 
dling the black hole which forms in the collapse of a neu- 
tron star fl5l 20l 2 111. a nd in binary systems in general (see 
e.g. iTTeL frTfmlH^l and references therein). In d, in 
particular, it was first shown with three-dimensional (3D) sim- 
ulations that when using the puncture gauges, no special treat- 
ment beyond standard artificial dissipation for the metric vari- 
ables is necessary to follow stably the black-hole formation 
and evolution. 

In this paper we show that the numerical evolution, with 
the puncture gauge, of a collapsing spherical star approaches 
the trumpet solution at late times. This result is not trivial for 
a number of reasons. First, because in the collapsing space- 
time, and in particular in the portion filled by matter, there is 
no time-like Killing vector and it is not clear how a stationary 
end-state can be found. Second, at the continuum level, the 
two spacetimes are clearly different since, for example, slices 



in the two foliations have a different topology. Even restrict- 
ing our attention to the vacuum region of the collapsing space- 
time and assuming a Killing slice compatible with the "1+log" 
lapse condition, it is not obvious that the resulting slice will be 
exactly that reported in u3f] . Finally, the collapsing matter in- 
side the apparent horizon is observed to disappear but it is not 
clear whether this reflects a physical or a numerical behavior. 

We therefore seek answers to the following questions: 
(i) How well do the two endstates agree? (ii) Where does 
the matter go after the apparent horizon is formed? (iii) How 
do the gauges find a stationary slicing of a spacetime without 
an explicit time-like Killing vector? In order to address these 
questions we evolve the spacetime of a star collapsing to a 
black hole using a "1+log" slicing condition and two variants 
of the Gamma-driver shift condition. The different results are 
then compared quantitatively with the trumpet solution, and 
spacetime diagrams are constructed to follow the motion of 
the matter. Overall we find that the two spacetimes tend to a 
common stationary solution at late times and that this agree- 
ment is primarily caused by the shift condition, which pushes 
grid-points away from regions of high curvature, preventing 
the matter from being resolved on the numerical mesh. 

The plan of the paper is as follows. In Sect. [II] we briefly 
review our numerical methods, presenting the numerical re- 
sults and their interpretation in Sect. [Ill] Our conclusions are 
summarized in Sect. [IV] while we dedicate the Appendix [A] 
to a description of our method for constructing spacetime di- 
agrams for spherical spacetimes. 



II. NUMERICAL SETUP 

We perform numerical simulations both in explicit spher- 
ical symmetry by means of the one-dimensional (ID) code 
described in J3 and in 3D by means of BAMMATTER JH]. 
The latter is developed extending the h ydro dynamics solver of 
the spherical code into the BAM code 12611 . More specifically, 
we solve the full set of Einstein equations in the 3 + 1 formal- 
ism coupled to general relativistic hydrodynamics (GRHD). 
In our spherically-symmetric simulations we adopt both the 
BSSNOK and Z4c formulations of the Einstein field equations 
(see JH for details), while in 3D we use only the BSSNOK 
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formulation. 

As mentioned above, our gauges are built upon the "1+log" 
condition (jl for the lapse a and the Gamma-driver condition 
for the shift /3* H written in the form 



d t a 



(1) 

(2) 



where we always choose Hl =2/ a, with either [is = 1 or 
ps = a 2 and rj = 2/M. The Gamma-driver condition (f2]i 
may be obtained as the first integral of what is implemented in 
most numerical-relativity codes. A comparison of the behav- 
ior of the condition © with the standard form in the evolution 
of puncture black holes is presented in |4j]. 

For the evolution of the matter we adopt the Valencia flux- 
conservative formulation of the equations of GRHD I27U29I1 
for a perfect fluid. A property of flux-conservative formu- 
lations is that they preserve (exactly in the continuum) cer- 
tain integral quantities such as the total rest-mass of the fluid, 
Mo, or its momentum and energy density. The equations of 
state (EoS) employed to describe the fluid are the ideal-gas or 
the poly tropic EoS l28ll : no significant differences were found 
when using one or the other. 

In both codes the numerical evolution in time relies 
on the method-of-lines with Runge-Kutta integrators and 
finite difference approximations. The GRHD equations 
are solved by means of a high-resolution shock-capturing 
(HRSC) scheme based on the Local Lax-Friedrich (LLF) 
central scheme l30ll combined with Convex-Essentially-Non- 
Oscillatory (CENCOjp] reconstruction JMS. 

In analogy with IU5I1 and in contrast to II2411 . we do not find 
it necessary to excise the hydrodynamical variables in the col- 
lapse evolutions. Artificial Kreiss-Oliger dissipation is added 
in the evolution of the metric fields following the standard pro- 
cedure i2oTl . Instead, we find it important to set the GRHD 
eigenvalues to zero if unphysical values are computed ("eigen- 
value excision"). Note that this is compatible with the use 
of the LLF central scheme, which requires only an estimate 
of the local speed. Unphysical values can be produced in a 
neighborhood of the center of the collapse, due to numerical 
errors. In 3D simulations we also found it important to set a 
ceiling on the Lorentz factor W in order to prevent the code 
from crashing after the formation of the apparent horizon. In 
particular W is set to the ceiling value W ce n = 10 10 when 
(and only when) the velocity becomes larger than the speed 
of light. Simulations employing the excision of the hydrody- 
namical quantities were also performed. In this case, the mat- 
ter variables are set to the atmosphere value in a small region 
well inside the apparent horizon. As expected, no differences 
were found. In the following we focus on simulations without 
excision. 

The initial stellar model is an unstable spherical configura- 
tion widely used in the literature, e.g., [28] . Adopting units 
in which c = G = Mq = 1, the configuration chosen has 
central rest-mass density p c = 7.993 x 10~ 3 , gravitational 
(ADM) mass M = 1.448 (M = 1.535) and circumferential 
radius R = 5.838 (isotropic coordinate radius r# = 4.268). 
The collapse is triggered imposing a negative velocity per- 
turbation which is larger than the truncation error rf28tl . For 
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FIG. 1: Radial profiles of the rest-mass density p (top panel), of the 
lapse function a (middle panel), and of the shift /3 = /3 r (bottom 
panel) at some representative times of the evolution. Note that at 
t ~ 50 M an apparent horizon is first found and that by t — 300 M 
the amount of matter on the final time slice is essentially that in the 
atmosphere. 



the vacuum simulations we use standard puncture data II 1311 
in isotropic coordinates and with an ADM mass which is the 
same as that of the star. 

The simulations in ID were performed on a grid with uni- 
form spacing with resolutions A r = 0.02,0.01,0.005. The 
simulations in 3D, on the other hand, were performed impos- 
ing an octant symmetry on a cell-centered Cartesian grid with 
8 fixed mesh refinement levels, with the resolutions of finest 
level given by A xyz = 0.05, 0.03125, 0.025, and where the 
resolution doubles from one level to the next. 



III. RESULTS 

We next describe the results of the different numerical sim- 
ulations. Before that, however, it is useful to recall some of 
the basic properties of the evolutions that from the puncture 
spacetime lead to the trumpet solution or that from a stellar 
spacetime lead to a black-hole solution. 

Puncture spacetime. The evolution of the puncture space- 
time is described in Ref. ll3l [Till . The initial Schwarzschild 
slice in isotropic coordinate evolves, driven by the puncture 
gauges, to the stationary trumpet solution. Approximate sta- 
tionarity is achieved within an evolution time of about 60M, 
depending on the details of the gauge choice. Asymptoti- 
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cally, there is a coordinate singularity at the puncture point 
at r = 0, which corresponds to a sphere with Schwarzschild 
radius Rq w 1.3M. The trumpet slice extends from Rq to 
spatial infinity, i°, while the wormhole slice of the initial data 
reaches from the outer i° to an inner spatial infinity. Con- 
ceptually it is important to distinguish between the analytical 
solution and the numerical solution on a discrete grid, and also 
between the wormhole slice of the analytical evolution and the 
trumpet slice which is only reached asymptotically. The evo- 
lution starts with a wormhole slice which, due to continuity of 
the analytical solution, remains a wormhole slice. It asymp- 
totes to a trumpet slice in the sense that the region inside i?o 
has a coordinate size that tends towards zero. Numerically, 
this inner region is effectively excised once its coordinate size 
drops below the grid spacing. Assuming that there is no grid 
point at the puncture itself, after a short evolution time the 
numerical grid only has grid points outside i?o, that is only 
the trumpet part of the initial wormhole is represented. An- 
ticipating the discussion of matter, the key difference is that 
the initial data with matter lives neither on a wormhole nor a 
trumpet slice, since the matter "covers" the inner region of the 
slice. 

Collapsing spacetime. The dynamics of a collapsing un- 
stable star in general relativity has been discussed in great de- 
tail in a number of p apers and we refer the interested reader 
to Refs. 01 HI 29] for some of the most recent work. For 
the test-case considered here (jug = 1), the introduction of 
the perturbation is sufficient to drive the star over the stabil- 
ity threshold and induce its gravitational collapse. As a re- 
sult, matter essentially freely falls towards the center, leading 
to an exponential growth of the rest-mass density and a re- 
lated rapid variation of the metric functions. This is shown 
in the top panel of Fig. Q] which exhibits the radial profile 
of the rest-mass density at some representative times. The 
"singularity-avoiding" properties of the "1+log" slicing con- 
dition drive the lapse function to very small values near the 
origin. This is shown in the middle panel of Fig.Q] while the 
bottom one refers to the evolution of the shift. Note that as 
the matter rushes towards the center, the shift still succeeds in 
arresting the motion of the radial coordinates outside the mat- 
ter distribution (ignoring the atmosphere), thereby preventing 
slice stretching (in the sense of stretching of the spatial coor- 
dinates) as for the black hole evolutions. At about t ~ 50 M, 
an apparent horizon forms indicating unambiguously the pres- 
ence of a black hole. This time represents the time when a 
first comparison between the two spacetimes can in principle 
be made. We also note that because the coordinate radius of 
the apparent horizon is initially r^n < 2 M, part of the matter 
is outside of it, but it is then rapidly accreted. 

As highlighted in 11511 . the gauge conditions in Eqs. (Q~|i-(|2]i 
allow us to follow the subsequent evolution without having 
to excise either the hydrodynamical or the gravitational field 
variable. The conservation of the rest-mass is very good up 
to the horizon formation: AM o (t)/M o (0) < 0.05%. Af- 
ter the formation of the apparent horizon, the matter inside 
is observed to disappear from the numerical grid, so that by 
t = 300 M the amount of matter on the slice is that in the 
atmosphere, Mq(300) ~ 10 -6 Mq(0). As we will discuss in 



Sect. 1111 Bl a suitable change in the value of /xg for the shift 
condition can avoid this behavior. 

As a side remark we note that in our ID simulations we 
observe a better behavior of the Z4c formulation of the Ein- 
stein equation with respect to the BSSNOK one in terms of 
constraint violation and long-term stability. In particular, at 
t = 2500 M, the L2 norm of the Hamiltonian constraint is 
about 4 orders of magnitude lower with Z4c. Furthermore, at 
t = 2500 M, the irreducible mass of the final black hole in the 
Z4c evolutions is within the numerical error of the ADM mass 
of the initial data. In contrast, the BSSNOK simulations dis- 
play a significant deviation in the irreducible mass of the black 
hole after t = 300 M, In view of this, in what follows the ID 
results we will present refer exclusively to those obtained for 
the Z4c formulation. 



A. Agreement of the end-states 

In this Section we investigate the asymptotic slice reached 
at the end of the evolution by the collapsing star and com- 
pare it with the corresponding trumpet one. This is shown in 
Fig. 12 which reports several metric fields at time t = 300 M. 
This time is well after the apparent horizon is first found 
(i.e., tAH ~ 50 M), and represents in both cases a time when 
the solution has become essentially stationary. From top to 
bottom the different panels refer respectively to: the trace of 
extrinsic curvature (K), the conformal factor (%), the lapse (a) 
and the shift (/3 r ). It is quite apparent that at the selected time 
the two spacetimes are extremely similar and it is difficult to 
distinguish the solutions two by a visual inspection. 

To obtain a more quantitative estimate of the differences 
we have computed the behavior of the fields in the collapsing 
spacetime near the origin and obtained the following fitting 
functions for the stationary solution for K< 1: 
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The fits for (01 and © contain the exponent of r as a fitting pa- 
rameter, but not for Q. The result agrees well with the corre- 
sponding expressions in rf32"tl . In particular, the non-integer ex- 
ponent for the lapse in (O is very close to the analytic expres- 
sion for the trumpet solution in isotropic coordinates, which 
has exponent 1.091 ll32ll . This is interesting since the numer- 
ical coordinates are not isotropic, although one could argue 
that close to the puncture the isotropy of the initial data is 
maintained during the evolution. 

This result is confirmed when considered also in a 
coordinate-independent manner. Following the prescription 
suggested in |3j], we analyze the dependence of the lapse ver- 
sus the extrinsic curvature and report in Fig.[3]the differences 
at different times. More specifically, we show with solid lines 
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FIG. 3: Relative difference Aa/a = a punc /a co n — 1, between the 
lapse of the puncture evolution, a punc , and that of the collapsing 
star, a co n, as a function of the trace of the extrinsic curvature, K, at 
different coordinate times. For a fair comparison, data at the same 
K are found by interpolation. The dashed lines refer to the relative 
differences between the analytic result and the puncture evolution. 
By construction, for early times this difference is much smaller than 
the difference between puncture and matter simulations. 
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where R is the Schwarzschild radius, a' is the derivative of 
a with respect to R and we set M = 1. As expected, the 
relative difference is in this case much smaller initially, but it 
becomes comparable with the one computed for the collapsing 
spacetime at later times. 

It is possible to model the behavior of the data from the 
collapsing spacetime near the origin and at t = 300 M as 



FIG. 2: Radial profiles of representative field variables at time t = 
300 M from the ID simulations of the puncture and the collapsing 
star. From top to bottom: trace of extrinsic curvature (K), conformal 
factor (x), lapse (a) and shift (/? = /3 r ). 



lapse of the puncture evolution, a pun c, and that of the col- 
lapsing star, a co ii, either when the apparent horizon has just 
formed (t ~ 50 M) or when the solutions have reached a sta- 
tionary stage (t — 300 M). 

It is clear that the relative difference decreases in time and, 
at time t = 300 M, it is below 0.1%. By performing conver- 
gence tests we have also determined that the numerical errors 
are at the same level as the 0.1% disagreement. Also reported 
in Fig.[3]with dashed lines is the relative difference between 
the puncture data and the analytic solution for the trumpet so- 



KM~ 0.30 -0.92 a. (7) 

Similarly, a Taylor expansion of Eq. © around a = (i.e., 
for the values of the lapse near the puncture) can be performed 
by using the implicit function R(a) in l32il . yielding 

K(a) = 0.300937- 0.930916a + O (a 2 ) , (8) 

which closely agrees with Eq. (0. 

Finally, we note that we did not find significant differences 
between the ID and the 3D results for the collapsing space- 
time. This is summarized in Fig.|4j where the top panel shows 
a versus K for the ID and the 3D data at time t = 300 M: 
they are visually indistinguishable. Similarly, the bottom 
panel shows that the relative differences between the ID and 
the 3D data are generically below 2%. 
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FIG. 4: Slicing from ID and 3D simulations. Top panel: the lapse 
(a) is plotted versus the trace of the extrinsic curvature K at time 
t = 300 M for the collapse in ID (red solid line) and the punc- 
ture evolution (blue solid line) spacetimes, and the collapse in 3D 
(red dashed line) and collapse (blue dashed line) spacetimes. Bottom 
panel: relative differences. Note that ID results, computed on a finer 
grid, are interpolated for the comparison. 
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FIG. 5: Evolution of total rest-mass, Mo(t), (red dashed lines) and of 
the irreducible mass, Mi IT (t), (blue solid lines) normalized respec- 
tively to the initial value and to the ADM mass. The top panel refers 
to a simulation with gauge speed /i 3 = 1, and the bottom panel to 
evolutions with gauge speed fi s — a 2 . 



B. Where does the matter go once inside the horizon? 

As discussed above, the numerical evidence obtained when 
using the Gamma-driver shift condition with /ig = 1 is that 
the matter inside the apparent horizon is progressively "dis- 
sipated" (see also lfl5ll where this was first discussed). This 
behavior is clearly displayed in the upper panel of Fig. [5] in 
which the total rest mass normalized to the initial value is 
plotted in time (dashed red line). Note that at about 10M af- 



ter the formation of the apparent horizon the normalized rest 
mass drops smoothly to zero except for a small bump near 
t = 65Af generated by numerical errors inside the horizon. 
Also shown is the irreducible mass of the black hole normal- 
ized to the ADM mass (solid blue line), which is obviously 
zero before the apparent horizon is found. This behavior may 
appear puzzling since it is geometrically clear that the mat- 
ter cannot leave the foliation (see e.g., Fig. 6.11 of |[33tl ). In 
addition, since we are using HRSC schemes, our numerical 
methods should be sufficiently robust, for reasonable spatial 
resolutions, to handle extremely large gradients as the matter 
piles up in the collapse. Indeed, we have experimented with 
various reconstruction algorithms and resolutions and found 
that both affect slightly the rate and the initial time of the dis- 
appearance of the fluid. More specifically, more dissipative 
schemes result in an earlier disappearance of the matter, while 
higher resolutions can delay it. However, neither improve- 
ment prevents the disappearance of the matter from the grid. 

In order to establish whether this behavior is instead due 
to an excessive stretching of the spatial coordinates generated 
by the Gamma-driver shift condition, we have performed the 
same simulations using either a shift gauge speed fis = a 2 , 
or simply setting /?* = 0. We recall that this is possible since 
the "1+log" lapse condition is a pure slicing condition, so that 
the foliation is unaffected by a change of radial gauge (see for 
example HUH). 

The results of these tests are shown in the bottom panel 
of Fig. [5] where we report with the evolution of the total 
rest-mass (dashed red lines) and of the normalized irreducible 
mass (solid blue lines) when /is = a 2 . Clearly, while the 
behavior of the irreducible mass is independent of the choice 
for us, that of the rest-mass is not. When using the /us = 
radial gauge, in fact, the matter remains on the numerical grid, 
so that the rest-mass is conserved well beyond the formation 
of the apparent horizon. The cause of this difference must 
therefore to be attributed to the large stretching of the spatial 
coordinates with /j,g = 1. Specifically, we identify two ef- 
fects. Firstly matter falls inside the innermost grid point as 
the inner boundary is effectively an outflow boundary for the 
matter; Mo starts to drop to the atmosphere value approxi- 
mately 10 M before the stellar surface passes through the in- 
nermost gridpoint. At this time, the areal radius of the in- 
nermost gridpoint grows rapidly to approximately 1.9 M (see 
also Fig. [6] later). Secondly the effective resolution in areal 
radius near the time of the stellar surface passing through the 
innermost gridpoint is approximately an order of magnitude 
lower than in the initial data. The stretching is thus so large 
that the matter "percolates" through the grid as the numeri- 
cal methods are not able to reproduce its steep gradients. As 
a result, the spacetime "empties" itself and this explains the 
very good match with the trumpet solution reported in Fig. [3] 
Conversely, when /is — a 2 , the radial gauge does not distort 
the grid significantly, allowing for an excellent conservation 
of the rest-mass on the grid. Note that the latter is not the re- 
sult of the rather high spatial resolution, but it is simply the 
result of the coordinate time "freezing" induced by the col- 
lapsed lapse function. Of course, because the spacetime is not 
able to remove its matter content, the match with the trumpet 
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solution is in this case much worse and Act /a ~ 10 for 
KM > 0.21. 

Not surprisingly, much of what we discussed so far for 
[is = o? applies also when considering /3 l = 0. However, 
the same dynamics that leads a f3 l = shift condition to fail 
in a curved spacetime, is responsible for the late-time failure 
of the simulations having [is = o? as gauge speed. As a re- 
sult, a Gamma-driver shift condition with [is = 1 appears to 
be the most robust choice for all those situations in which a 
compact fluid object may collapse to a black hole. 

Finally we show in Figs. [6] and [7] the spacetime diagrams 
as built from the numerical data. The axes in Fig. prefer to 
the time and spatial coordinates used respectively in the sim- 
ulations with [is = 1 (upper panel) and [is = ce 2 (bot- 
tom panel). The shaded green area corresponds to the region 
of the spacetime covered by matter, vertical lines are line of 
constant Schwarzschild radius while horizontal lines are lines 
of constant coordinate time. The thick red line emerging af- 
ter t ~ 50 M is the apparent horizon. Fig. [6] can be com- 
pared with other numerically-generated diagrams of collaps- 
ing spacetimes obtained either in other gauges ll35tl . or with 
excision tec hniq ues 112911 . or for puncture evolutions of single 
black holes ll ill . The figure clearly shows how the matter is 
"squeezed" from the numerical grid when [is = 1, while it 



remains on the grid when [is 



In both cases an appar- 



ent horizon is found. The comparison of the lines of constant 
Schwarzschild radius in the two panels highlights the stretch- 
ing of the spatial coordinate discussed before. 

Figure Q displays the causal spacetime diagram of the col- 
lapse. Similarly to Fig. |6]the thick red line denotes the appar- 
ent horizon, the thick green line the surface of the collapsing 
star and the shaded area shows the region filled by the matter. 
The coordinate (X, T) are constructed so that the light speeds 
have magnitude one, with light cones opening at 45 degree 
in the diagram. The method employed for the construction 
is discussed in Appendix [A] Because easier to use, the data 
employed to construct the diagram refers to a ID simulation 
with [is = a 2 , but in principle we expect that a very similar 
diagram would be produced with [is = 1. At early times the 
causal spacetime diagram clearly shows the extremely high 
speed of the stellar surface and the causal separation between 
the points initially belonging to the exterior and those that are 
causally connected to the collapse (see the gauge-wave prop- 
agating at the speed of light). At late times it is evident that 
the collapse of the lapse prevents the slices from evolving for- 
ward into the singularity; furthermore the inset demonstrates 
that there is always some numerical data inside the horizon. 




FIG. 6: Spacetime diagram of the collapsing star. Data are from 
ID simulations with gauge speed (is = 1 (top panel) and [is = 
a 2 (bottom panel). The horizontal blue lines are lines of constant 
coordinate time. The thick red line denotes the apparent horizon. 
The vertical blue lines are lines of constant Schwarzschild radius R 
which values are on top of the lines. The shaded green area bounded 
by the thick green lines shows the region of the matter. 



IV. CONCLUSION 

By comparing the numerical evolution of a single puncture 
with that of a collapsing star when using the same puncture 
gauges, we have shown numerically that the two spacetimes 
tend to the same trumpet solution at late times, a possibility 
already conjectured in 11511 . Let us explicitly address the ques- 
tions raised in the introduction: (i) In the domain covered by 
the numerical coordinates, at late times, the spacetimes agree 



within the precision of our calculations, (ii) The apparently 
bizarre agreement, which is not possible at the continuum 
level, is caused primarily by the Gamma-driver shift condi- 
tion, which stretches the numerical spatial grid. The matter 
is rapidly forced inside the innermost grid point, preventing 
it from being resolved and effectively removing it from the 
spatial slice. Thus in the domain covered by the numerical 
coordinates the spacetimes agree, solving the apparent con- 
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are given by 



FIG. 7: Causal spacetime diagram of the collapsing star. Data are 
from ID simulations with gauge speed fis = ■ The coordinate 
axes are constructed so that the light speeds have magnitude one. 
The horizontal blue lines are lines of constant coordinate time. The 
thick red line denotes the apparent horizon. The shaded green area 
bounded by the thick green lines shows the region of the matter. 



tradiction. (iii) Since the matter is lost inside the innermost 
gridpoint, at late times, on the numerical grid the evolution is 
simply that of the vacuum Schwarzschild solution. 

We also note that if the matter in the slice is sufficiently 
compact and in the vacuum region and one has a Killing slic- 
ing which is compatible with the "1+log" gauge, then in the 
vacuum region the slice must agree with the original "1+log" 
trumpet Ool . Our analysis also demonstrates that, whilst the 
puncture gauges with fig = 1 results in robust numerical evo- 
lutions in a collapse scenario, by construction the resulting 
coordinates are not appropriate for a detailed study of the dy- 
namics of the matter near the singularity. 



Appendix A: The construction of spacetime diagrams 

In this appendix we describe our approach to the construc- 
tion of spacetime diagrams in spherical symmetry. The ap- 
proach is based on that of l37ll . The most general spherical 
line element can be written as 



ds 2 = -a 2 dt 2 + lrr {Ar + f3 r dtf + 7ee dft 2 , 



(Al) 



where dO 2 = d9 2 + sin 2 9 d<p 2 is the standard metric on the 
two-sphere. The incoming and outgoing radial light speeds 



-/3 r ± 



a 



(A2) 



Ingoing and outgoing null coordinates (u, v) satisfy the equa- 
tions of motion 



d t u 



-c+d r u, dtv = —C—d r V, 



(A3) 



in coordinates (t,r,6,(f)). Assuming that near the outer 
boundary r* space is almost flat, we arrive at the boundary 
condition 



v(t, r*) = 2t — u. 
We introduce the scalar fields 



X = \(v-u), 



T = -(u + v). 



(A4) 



(A5) 



They satisfy the equations of motion 



d t X = f3 r d r X 



-_d r T, d t T = p r d r T 



-_d r X. 
(A6) 



and are naturally adjusted to the causal structure of the space- 
time. In terms of (T, X, 9, 4>) the spherical line-element be- 
comes 



ds 2 



^(-dT 2 + dX 2 ) + "fggdn 2 



where i\> is a conformal factor determined by 

/ 4 _ Irr 



{d r xy ~ {d r Tf 



(A7) 



(A8) 



To construct spacetime diagrams we simply evolve the fields 
(u, v). As initial data we choose u = —r and v = r. There 
are potential problems associated with the breakdown of the 
coordinate chart, but we will not concern ourselves with those 
issues here. From our computations we observed however that 
the method requires further development (e.g., a treatment of 
the fields at the horizon) to construct spacetime diagrams for 
every numerically generated spacetime (cf. Sec. Ill B of [ 1 1 ]). 
In particular we found the construction of causal spacetime 
diagrams with fig = 1 troublesome, perhaps because more 
sophisticated boundary conditions are required. 
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